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Ch . Abstract 

The results of numerical experiments on chaotic ('turbulent') dynamics of 
^ ■ the second sound in helium II are presented and discussed based on a very 

simple model proposed and theoretically studied recently by Khalatnikov and 
Kroyter. Using a powerful present-day techniques for the studying nonlinear 
phenomena, we confirm their results on the stationary oscillation in helium 
U: and its stability as well as on a qualitative picture of successive transitions to 

C I limit cycles and chaos. However, the experiments revealed also a much more 

complicated structure of the bifurcations than it was expected. The fractal 
structure of chaotic attractors was also studied including their noninteger 
dimension. Surprisingly, a very simple model used in all these studies not 
^ ' only qualitatively represents the behavior of helium in laboratory experiments 

OO . but also allows for a correct order-of-magnitude estimate of the critical heat 

I pumping into helium at bifurcations. 

o 

g ; 1 Introduction 

The present studies were stimulated by the theoretical analysis of the dynamics 
. of the second sound in helium driven by the external heat periodic perturbation 

J> ! [HI]. In turn, this theory was developed to explain the experiments on the second 

I sound propagation in a liquid helium with a free surface 0, Originally, such 

}-j I experiments were intended for a new scheme of holography using the deformation of 

the helium surface under the standing waves in the bulk of the liquid 0]. To the best 
of our knowledge, such a project was never realized. The reason is very simple: a 
negligible surface deformation which is determined by the oscillation of the pressure 
in a second-sound wave, which is much smaller than the temperature oscillation. 
Nevertheless, the experiments were continued, and an interesting dynamics was 
observed, namely the transition, at a sufficiently high heat pumping into the helium, 
from a stationary second-sound wave to its chaotic motion, or turbulence 0. This 
was just the subject of the theoretical analysis in Ref.|^. The latter was necessarily 
restricted, even in the simplest model considered, to the calculation of the stability 
of a stationary standing wave. For this reason the authors [0 began also some 
numerical experiments which qualitatively confirmed their conjecture on a transition 
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to limit cycles and chaos. The main purpose of our numerical studies is to continue 
such experiments for a more detailed and reliable analysis of those transitions. 



2 The model 

We adopted the model developed in Ref.|jl| which can be specified by the following 
effective non-Hermitian Hamiltonian (in notations of Ref . with only a few minor 
modifications): 

H{ai, a2, t) = {ui - i7i)|air + (^2 - «72)|a2p + {Xafal + fale'''^^ + c.c.) (2.1) 

This Hamiltonian describes the two wave modes (oscillators) via complex phase- 
space variables aj [j = 1,2) and the complex small-oscillation frequencies Uj — 
ijj with phenomenological dissipation parameters jj. Two other parameters of 
the model, representing some effective nonlinear coupling of the two modes A and 
the driving 'force' (external heat pumping) / of frequency u, can be chosen real 
and positive. The model parameter / is related to the input heat power by an 
approximate expression 

dH 

P ^ = -2/^Im(ai) (2.2) 

provided the parameters 7^- and / are sufficiently small. Here aj are the new variables 
('slow' amplitudes of the wave modes) determined by the canonical transformation 

which eliminates the explicit dependence on time in Eq.(2.1). Below we omit the 
upper bar in all aj. Notice that the new Hamiltonian still exactly describes the 
original model (2.1). The corresponding motion equations have the form: 



icii = (Ai — 271)01 + 2Aa^a2 + / 
id2 = (A2 - 272)02 + Aa^ 



(2.4) 



where Ai = ui — u and A2 = 0^2 — 2a; are the detunings of the small-amplitude 
oscillation frequencies with respect to that of the driving force. 

The described model is actually the simplest possible one for representing a 
complicated dynamics of the real physical system under consideration. Particularly, 
only the lowest (third-order in Hamiltonian) nonlinear term is included. Meanwhile, 
the next, fourth-order term is known to give a comparable contribution, particularly 
to the nonlinear shift of the oscillation frequencies. In Ref. |1]] an argument is given 
for neglecting this term. In our opinion, the main argument is in that the next 
term does not lead to new qualitative effects, and moreover retains the order-of- 
magnitude of the description which is the best one could expect from such a primitive 
model. In any event, following Ref.|jll| we also restrict our studies to the equations 
(2.4). 
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3 Numerics 



In numerical experiments the set of equations (2.4) was integrated by itself (for 
trajectories) or together with the linearized equations (for Lyapunov exponents) 
using the Fehlberg version of the forth-order Runge - Kutta algorithm. A typical 
integration step was At ~ 2^^*^. The accuracy of integration was checked, as usual, 
by a change of the step size. In case of a stable motion (stationary oscillation, limit 
cycle) the criterion of accuracy was a negligible change of the trajectory. However, 
for chaotic motion this was never possible, and we checked instead the accuracy of 
statistical properties of the motion, for example of the Lyapunov exponents. 

We computed all (4) exponents A„ (Amax = Ai > A2 > A3 > A4) using also the 
standard method (see, e.g., Ref.|]5|). One of them, whose eigenvector goes along the 
trajectory, is always zero while the sum of all 

4 

Ar = ^ A, = -2(71 + 72) = const (3.1) 

n=l 

is the rate of the phase space volume contraction. Both conditions were also used 
to check the integration accuracy. 

All the model parameters but / were fixed to (in CGS units) 

Ai = A2 = -1500 
71 = 30 72 = 120 
A = 5400 and 
uj = 8000 

as in Ref.[jl|. The principal parameter A was calculated in Ref. 0. The rest were 
chosen close to those in a typical laboratory experiment ^. We add also the driv- 
ing frequency uj P], absent in Eq.(2.4), which we will need below for a quantitative 
comparison with the laboratory experiments. The remaining parameter /, related 
to the input heat power (see Eq.(2.2)), was the only varying one in this first series 
of our numerical experiments. 

In agreement with laboratory experiments the initial conditions were fixed at 

ai(0) = 02(0) = (3.3) 

which is important for the interpretation of numerical experiments below. 

The main results of our experiments are presented in Fig.l as a dependence of 
the maximal Lyapunov exponent Amaxif ) on the driving parameter /. 

Three qualitatively different regimes of the motion are clearly seen: 

(i) Amax < and, hence, all A's are negative. This means a stationary standing 
wave according to the original Hamiltonian (2.1) or the constant amplitudes 
{ttj =const) in the reduced motion equations (2.4) that is a fixed point in the 
amplitude space. This is the simplest regime well recognized in both laboratory 
and numerical experiments. 
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(ii) Krnax = 0. This may be a periodic (limit cycle) or quasiperiodic attractor. 
In the former (simplest) case the trajectory is closed. One way to recognize 
the cycle is simply to see a picture of the trajectory, or rather of some two- 
dimensional representation of the four-dimensional trajectory. This was done 
in Ref . |1[ , and we will not repeat it here. However, in case of a quasiperiodic 
attractor this visual analysis doesn't work. 

(iii) Amax > 0. This is the most simple and reliable criterion for chaotic motion. 
The visual difference between the chaotic attractor and a complicated limit 
cycle is not always clear and may be deceptive. 




Figure 1: The global dependence of the maximal Lyapunov exponent Amax on the 
driving force / revealing the three types of attractor in model (2.4) with parameters 
(3.2) and initial conditions (3.3). Empirical points are connected by line to guide 
the eye. 



As the force / grows the first very sharp bifurcation (fixed point/limit cycle) 
occurs at /i ^ 97 but it does not mean the fixed-point attractor disappears for 
larger /. Indeed, both the theory [|l| and our numerical experiments show that the 
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Figure 2: Same as in Fig.l for a very narrow interval of /, indicated by arrow, where 
a deep drop of A^aa; > occurs. 



attractor retains up to / ~ 380 for some initial conditions of the motion which are 
generally different from (3.3). What actually happens in the first bifurcation is the 
escape of the initial point (3.3) out of the basin of attraction Particularly, this 
explains the 'window' with a fixed-point attractor at larger / ^ 137 — 185 (Fig.l). 

At still greater / = /2 ~ 461 a new and much more interesting bifurcation 
occurs (limit cycle/chaotic attractor). Actually, it is a complicated chain of many 
bifurcations between periodic and chaotic attractors. This is especially clear in Fig. 2 
where a narrow part of dependence Amaxif) is shown in detail (cf. Fig.l). 

The structure of each chaotic attractor is rather intricate, and particularly frac- 
tal. An important characteristic of such a structure is the dimension which is gen- 
erally noninteger (see, e.g., Ref. 0). The dimension can be calculated from the 
Lyapunov exponents using the relation 

d = m+ ^"=^ 7 (3.4) 

|^m+l| 

where m is the largest integer for which the sum Sm = J2^=i-An > 0. This is the 
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simplest and widely used method for calculating the dimension. For example, on 
a fixed-point attractor (all A's are negative) m = 0, and d = 0. Similarly, for a 
periodic attractor (limit cycle) m = 1, and d = 1 {a closed trajectory). In case of 
a quasiperiodic attractor m = = 2 or 3, so that the trajectory is asymptotically 
restricted to a phase-space surface of dimension 2 or 3, respectively. We tried but 
couldn't find any quasiperiodic attractor in the model under cosideration. 

Unlike this, the dimension of a chaotic attractor is fractal because two of the four 
Lyapunov exponents are independent. In principle, three different cases are possible: 
(i) A2 = (Ai > 0); (ii) A3 = (Ai > A2 > 0); and (iii) A2 = A3 = (Ai > 0). 
In our experiments we have always seen the first case only and never the two latter 
ones. 

In the experiments presented in Fig.l for / < 600 (Ai < 140 the attractor 
dimension was in the interval {2 < d < 3) corresponding to m = 2 and = Ai 
(see Eq.(3.4)). For larger / the dimension grows due to increase of m {Sm = S3 = 
Ai + A3 > 0), and can be represented in the form: 

,(S,) = 3-| = 3+^-^<4 (3.5) 

In agreement with the physical meaning of the dimension it never exceeds the upper 
bound dmax = 4. Whether d is actually approaching this bound is an interesting 
question. However, it makes sense for the model only, since at such a big / the model 
becomes certainly inadequate for the real physical system under consideration. 

The above analysis of attractors in the model was confirmed by the computation 
of the motion spectra and correlation functions. 



4 Discussion 



How strange it may seem the very simple Khalatnikov - Kroyter model considered 
above does describe some qualitative results of the laboratory experiments P, |^. 
The most important of those is the observed transition from a stationary standing 
wave to something which looks like a chaotic motion or wave turbulence. 

In numerical experiments three types of attractors are clearly seen (Fig.l): the 
stationary standing wave, limit cycles, and chaotic attractors. In laboratory exper- 
iments it is generally rather difficult to be sure that a particular attractor is truly 
chaotic one rather than a very complicated limit cycle, for example. Nevertheless, 
the experimental observations presented in Ref. 0, ^ provide some plausible indica- 
tions of the transition to turbulence. The most clear evidence of that can be found 
in Fig. 2 [Q] where the motion of the helium free surface is shown for a few (3) val- 
ues of the input heat power. According to these data, the first transition (stable 
surface/limit cycle) occurs at a power Pi 18 mW. The inequality is explained by 
the absence of experimental data for lower power, so it remains unclear whether the 
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transition at Pi ~ 18 mW is the first escape of the system out of the basin of attrac- 
tion or it is a window at a higher power hke in our Fig.l for the model. The second 
(most interesting) transition to an apparent turbulence takes place somewhere in 
the interval P2 ~ 18 — 93 mW. Here we take the lower P2 equal to the upper Fi 
for the same reason, the absence of intermediate experimental data. 

In the model, the critical power input can be estimated from Eq.(2.2) as 

Pi f5i 2.6 mW at f ^ 95, Im(ai) ^ -0.016 (4.1a) 

for the first transition and the first escape out of the attraction basin (see Fig.l) or 

Pi ^ 4.8 mW at / 170, Im(ai) -0.018 (4.16) 

for the second escape. As to the transition to a chaotic attractor, the critical power 
in the model is 

P2 fa 26 mW at / fa 450, Im(ai) fa -0.036 (4.1c) 

In such a simple model it seems to be a quite reasonable agreement (especially for 
the transition to turbulence), and in view of a considerable uncertainty in laboratory 
experiments. 

The success of the Khalatnikov - Kroyter model can be explained, at least partly, 
by the special choice of the parameters of experiments when one of the wave modes 
was in resonance with the input heat oscillation. Since the dissipation in higher 
modes is rapidly increasing, only a few lower modes are excited initially which turned 
out to be already sufficient for the transition to turbulence. For a higher power input 
the number of excited modes also rapidly grows, and the model becomes inadequate. 

An interesting question about the maximal dimension of a chaotic attractor in 
the model (see Section 3), being irrelevant for the physical model under consideration 
here, may still be of importance for the theory of dissipative dynamical systems, and 
deserves further studies. 
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